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We calculate the Schiff moment of the nucleus 199 Hg, created by 7rNN vertices that are odd under 
parity (P) and time-reversal (T). Our approach, formulated in diagrammatic perturbation theory 
with important core-polarization diagrams summed to all orders, gives a close approximation to 
the expectation value of the Schiff operator in the odd-A Hartree-Fock-Bogoliubov ground state 
generated by a Skyrme interaction and a weak P- and T-odd pion-exchange potential. To assess 
the uncertainty in the results, we carry out the calculation with several Skyrme interactions, the 
quality of which we test by checking predictions for the isoscalar-El strength distribution in 208 Pb, 
and estimate most of the important diagrams we omit. 



I. INTRODUCTION 



The existence of a permanent 1 electric dipole moment (EDM) in lcptons, neutrons or neutral atoms is direct 
£S| ' evidence for time-reversal (T) violation. Because of the CPT theorem, the search for EDMs can provide us valuable 
information about sources of CP violation. Though a phase in the Cabibbo-Kobayashi-Maskawa matrix is enough 
to account the level of CP violation in kaon and B-meson decays, it cannot explain the observed matter/ anti-matter 
asymmetry in the Universe. Physics that can (and new physics at the weak scale more generically) should produce 
EDMs not far from current upper limits. 
' So far no EDMs have been observed, but experiments are continually improving. Here we are interested in the 
, conclusions that can be drawn from the measured upper limit [| in 199 Hg, a diamagnetic atom. The largest part of 
its EDM most likely comes from T violation in the nucleus, caused by a T-violating (and parity-violating) component 
<^> of the nucleon-nucleon interaction. The atomic EDM is generated by the subsequent interaction of the nucleus with 
*0 ' the electrons. 

That interaction is more subtle than one might think. If the nucleus and electrons were non-relativistic point- 
particles interacting solely via electrostatic forces, the electrons would rearrange in response to a nuclear EDM to 
cancel it essentially exactly. Fortunately, as was shown by Schiff Q , the finite size of the nucleus leads to a residual 
atomic EDM. It turns out, however, that the relevant nuclear quantity is not the nuclear EDM but rather the nuclear 
"Schiff moment" 



S = (y \S z \*o) , (1) 



which is the nuclear ground-state expectation value, in the substate l^o) with angular momentum projection Mj 
; I ■ equal to the angular momentum </, of the z-component of the "Schiff operator" 



8 = To £ fa - I r P ■ (2) 

p— l ^ ' 

Here e is the charge of the proton, (r 2 ) c h is the mean squared radius of the nuclear charge distribution, and the sum 
is restricted to protons. 

For the Schiff moment to exist, P and T must be violated by the nuclear Hamiltonian. We assume that whatever its 
ultimate source, the T violation works its way into a meson-mediated P- and T-violating NN interaction generated 
from a Feynman graph containing a meson propagator, the usual strong meson-NN strong vertex and a (much 
weaker) P- and T-violating meson-NN vertex. The second vertex can take three different forms in isospin space. 
References 0, 0, Hi showed that short-range nuclear correlations and a fortuitous sign make the contribution of p- 
and w-exchange to the interaction small compared to that of pion-exchange if the T-violating coupling constants of 
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1 Permanent or static, rather than induced. 
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the different mesons are all about the same, and so we neglect everything but pion exchange. The most general P- 
and T-odd NN potential then has the form (fi = c = 1) 



2 - 

W(r a - r b ) = - J^- { [go (t„ ■ n) - f + t§) + g 2 (3r a V 6 * - r a ■ n )] (cr a - <r b ) 



87tton 

exp {-m^\r a - r b \ 



-f «-^)K + o- b )}-(r a -r 6 ) 



m w |r a - r 5 | 2 



1 



m„\r a - r b \ 



(3) 



where m,r is the mass of the pion, ton that of the nucleon, r z \p) = — \p), g = 13.5 is the strong 7rNN coupling constant, 
and the gj are the isoscalar (i = 0), isovector (i = 1), and isotensor (i = 2) PT-odd 7rNN coupling constants. A word 
of caution here: more than one sign convention for the g's is in use. Our go and gi are defined with a sign opposite 
to those used by Flambaum et. al [fJQ and by Dmitriev et. al HQ- 

The goal of this paper is to calculate the dependence of the Schiff moment of 199 Hg on the T-violating 7rNN 
couplings (we leave the dependence of these couplings on fundamental sources of CP violation to others) so that 
models of new physics can be quantitatively constrained. An accurate calculation is not easy because the Schiff 
moment depends on the interplay of the Schiff operator with complicated spin- and space-dependent correlations 
induced by the the two-body interaction W. In the early calculation by Flambaum, Khriplovich and Sushkov almost 
two decades ago @ , the correlations were taken to be admixtures of simple 1-particle 1-hole excitations into a Slater 
determinant produced by a one-body Wood-Saxon potential. 

More recent work 0, [ij made significant improvements by treating the correlations in the RPA after generating an 
approximately self-consistent one-body potential. However, that work used only the relatively schematic Landau- 
Migdal interaction (in addition to W) in the RPA and mean-field equations, and did not treat pairing self consistently. 
The reliance on a single strong interaction makes it difficult to analyze uncertainty. Such an analysis seems to be 
particularlyimportant in 199 Hg, the system with the best experimental limit on its EDM. The calculated Schiff moment 
of Refs. H,13 in that nucleus depends extremely weakly on the isoscalar coefficient go, a result of coincidentally precise 
cancellations among single-particle and collective excitations. They might be less precise when other interactions are 
used. 

Here we make several further improvements. Our mean field, which we calculate in 198 Hg before treating core 
polarization by the valence nucleon, includes pairing and is exactly self consistent. Pairing changes the RPA to 
the quasiparticle-RPA (QRPA), a continuum version of which we use to obtain ground-state correlations. Most 
importantly, we carry out the calculation with several sophisticated (though still phenomenological) Skyrme NN 
interactions, the appropriateness of which we explore by examining their ability to reproduce measured isoscalar-El 
strength (generated by the isoscalar component of the Schiff operator) in 208 Pb. The use and calibration of more 
than one such force allows us to get a handle on the uncertainty in our final results. 

The rest of this paper is organized as follows: Section II describes our approach and the Skyrme interactions we 
use, and includes their predictions for strength distributions that bear on the Schiff moment. Section III presents 
our results and an analysis of their uncertainty, including a calculation in the simpler nucleus 209 Pb that allows us to 
check the size of effects we omit in 199 Hg. Section IV is a conclusion. 



II. PROCEDURE FOR EVALUATING SCHIFF MOMENTS 



A. Method 



Our Schiff moment is a close approximation to the expectation value of the Schiff operator in the completely self- 
consistent one-quasiparticle ground state of 199 Hg, constructed from a two-body interaction that includes both a 
Skyrme potential and the P- and T-violating potential W. It is an approximation because we do not treat IF in a 
completely self consistent way, causing an error that we estimate to be small in the Section ITTT1 In addition, we do 
not actually carry out the mean-field calculation in 199 Hg itself. Instead, we start from the HF+BCS ground-state 
of the even-even nucleus 198 Hg and add a neutron in the 2p!/ 2 level. We then treat the core-polarizing effects of 
this neutron in the QRPA. A self-consistent core with QRPA core polarization is completely equivalent to a fully 
self-consistent odd-A calculation [ljj. We omit one part of the QRPA core polarization, again with an estimate 
showing its contribution to be insignificant. 

A good way to keep track of the two interactions and their effects is to formulate the calculation (and corrections 
to it) as a sum of Goldstone-like diagrams, following the shell-model effective-operator formalism presented, e.g., 
in Ref. The one difference between our diagrams and he usual "Brandow" kind is that our fermion lines will 

represent BCS quasiparticlcs rather than pure particles or holes. Our diagrams reduce to the familiar kind in the 
absence of pairing. 
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FIG. 1: First-order quasiparticle diagrams contributing to the Schiff moment. Diagrams (i), (ii), (iv) and (v) do not contribute 
if the valence nucleon is a neutron, which is the case in 199 Hg. Diagram (vi) is the complex conjugate of diagram (iii). 



We begin, following a spherical HF+BCS calculation in 198 Hg (in a 20-fm box with mixed volume and surface 
pairing fixed as in Ref. |l2(), by dividing the Hamiltonian into unperturbed and residual parts. The unperturbed 
part, expressed in the quasiparticle basis, is 

H a = T + Voo + Vu , (4) 



where T is the kinetic energy and V the Skyrmc interaction, with subscripts that refer to the numbers of quasiparticles 
the operator creates and destroys. The residual piece 2 is 



H [CS = W + V 22 + Vb + V31 + V 04 + V 40 



(5) 



The interaction W can also be expanded in terms of quasiparticle creation and annihilation operators; all the terms 
are included in H les , though Woo vanishes because W is a pscudoscalar operator. The valence "model space" of 
effective-operator theory is one-dimensional: a quasiparticle with u = and v — 1 (i.e. a particle, since it is not part 
of a pair) in the the a = (2p 1 / 2l m = 1/2) level. The unperturbed ground state |$ a ) is simply this one-quasiparticle 
state. Then the expectation value of S z , Eq. (QJ, in the full correlated ground state \^> a ) = \^q) is given by 



(*«|-S*|*„> = Af- 1 ^ 



l + H Tt 



Q 



Hr 



S Z 



Hn 



(0) 



Here e a is the single-quasiparticle energy of the valence nucleon, the operator Q projects onto all other singlc- 
quasiparticle states, TV is a normalization factor that, we will argue later, is very close to one, and the dots represent 
higher-order terms in H les . To evaluate the expression, we write S z in the quasiparticle basis as S z = Sn + 6*02 + S20 
(Soo vanishes for the same reason as Woo)- 

The zeroth-order contribution to the Schiff moment in Eq. © vanishes because the Schiff operator cannot connect 
two states with the same parity, and also because the Schiff operator acts only on protons while the valence particle 
in 199 Hg is a neutron. (There are no center-of-mass corrections to the effective charges ^3-) The terms that are first 



2 The BCS transformation makes V02 and V20 zero. 
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FIG. 2: Examples of diagrams contributing to the collective response to the Schiff operator. The sum of all these diagrams is 
represented in Fig. I^i). The broken line represents the action of the Schiff operator (as in Fig. 0, the zig-zag line represents 
the P- and T-violating interaction (also as in Fig. , and the looped line represents a generic Skyrme interaction. 



order in H les do not include the strong interaction V because it has a different parity from the Schiff operator. Thus 
the lowest order contribution to the Schiff moment is 



(Va\S Z \y a ) 



\ first— order 



w 



H 



s : - 



ill 



(7) 



where q\ is the creation operator for a quasiparticle in the valence level a and |— ) is the no-quasiparticle BCS vacuum 
describing the even-even core, so that |<I> Q } is just q\\— ). The contribution of Eq. J7J in an arbitrary nucleus can be 
represented as the sum of the diagrams in Fig. ^ the rules for which we give in the Appendix 3 . In 199 Hg, because 
the valence particle is a neutron, only diagrams (iii) and (vi) are nonzero. We can interpret diagram (iii) as the 
Schiff operator exciting the core to create a virtual thrcc-quasiparticlc state, which is then dc excited back to the 
ground state when the valence neutron interacts with the core through W . This diagram and its partner (vi) are 
what was evaluated by Flambaum et al. [||, though their mean field was a simple Wood-Saxon potential, their W 
was a zero-range approximation that didn't include exchange terms, and they neglected pairing. 

Core polarization, implemented through a version of the canonical-basis QRPA code reported in Ref. (with 
residual spurious center-of-mass motion removed following Ref. |l5j|'). can be represented by a subset of the higher- 
order diagrams. Because W is so weak, we need only include it in first order. The higher order terms in V that we 
include have the effect of replacing the two-quasiparticlc bubble in (iii) and (vi) of Fig. ^ with chains of such bubbles 
(see Fig. 0), as well adding diagrams in which the QRPA bubble chains are excited through a strong interaction of 
the core with the valence neutron. We therefore end up evaluating diagrams labeled A, Bl, and B2, in Fig. [3] (plus 
two more of type B in which the interaction W is below the Schiff operator). The explicit expression for diagram A, 
the first on the left in the figure, is 



(* a |S' 2: |* a )diag-A 



EEE 

A k>l k'>V 



15, 



)2 \kl}Z k , l , 



[ak'l'\W 3 i\a)£ 



-l 

A 



(8) 






(iii) diag-B2 



FIG. 3: QRPA diagrams contributing to the Schiff moment. The filled bubble represents an infinite sum of quasiparticle 
bubbles, including all the forward and backward amplitudes. The two B diagrams have partners (not shown) in which W acts 
below S z . 



3 These rules are similar to the ordinary Brandow/Goldstone-diagram rules but since fermion lines represent quasiparticles, their number 
need not be conserved at each interaction. In addition, we have only upward going lines because there are no quasiholes. 
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Here, Zh = Xb, + represents the QRPA amplitudes (the sum appears because the matrix elements of all our 
operators are real) and £\ is the energy of the collective state A in 198 Hg. The quasiparticle matrix elements (—\S§ 2 \kl) 
and (akl\W3i\a) are related to the usual particle matrix elements (k\S z \l) and (ak\W\al) through the transformations 
discussed in the Appendix. In the absence of QRPA correlations, the X and Y amplitudes are 1 or 0, and Eq. © 
reduces to that associated with diagrams (hi) and (vi) of Fig. ^ 
Diagrams Bl and B2 of Fig. have the explicit expressions 

(*a|5 Z |* a ) diag _Bl = -2^EE EH^l|^H(<I^W4'(*'iK|->(e.-e»)"V . (9) 

A i k>l k'yl' 

<*a|S 2 |* Q ) diag _B2 = ~ 2 EEE E^l^il J ) Z fe^-| 5 o 2 |fc0^(^'l^3i|a)(ea-e l )- 1 (fA-ea + e l )- 1 .(10) 

A i k>l k'>l' 

The factor 2 accounts for diagrams not shown in Fig. [3] in which W acts below the QRPA bubble, and the e's are 
quasiparticle energies. The difference between Eq. (JHJ) and Eq. (|10|) is mainly in the thrcc-quasiparticle intermediate 
states. 

A complete QRPA calculation that is first order in W would also include versions of diagram A in which W trades 
places with one of the Vs in the bubble sum. We don't evaluate such diagrams but estimate their size (which we find 
to be small) rom calculations in the simpler nucleus 209 Pb in the next section. We also use that nucleus to examine 
other low-order diagrams not included in the bubble sum of diagram A. 

Why do we expect the QRPA subset of diagrams to be sufficient? The reason is that they generally account for 
the collectivity of virtual excitations in a reliable way when calculated with Skyrme interactions. We illustrate this 
statement below with some calculations of isoscalar-El strength in 208 Pb. 

B. Interactions 

We carry out the calculation with 5 different Skyrme interactions: our preferred interaction SkO' |l2Lll6l (preferred 
for reasons discussed in Ref. 0])' an d the older, commonly used interactions SIII 0]j SkM* 0], SLv4 |20l|. and SkP 
|2ll | . To get some idea of how well they will work, we calculate the strength distribution of the isoscalar-El operator 

z N 

D o = E r l r p + E r « r « ■ ( n ) 

p— 1 n—l 

This operator is interesting because it is the isoscalar version of the Schiff operator (the isoscalar version of the second 
term in the Schiff operator acts only on the center of mass and so doesn't appear in Do). The isoscalar-El strength, 
measured, e.g., in 208 P b 1221 . s eems to fall mainly into two peaks. The high-energy peak, related to the compressibility 
coefficient K ^. I22L p3l 124 12511261 l2Sj . is observed to lie between 19 and 23 MeV, depending on the experimental 
method used |22,|23|. Recent interest has focused on a smaller but still substantial low-energy peak around 12 MeV, 
which has been studied theoretically in the RPA [2^, [2(| as well as experimentally |2^ ■ 

Figure 0] shows thepredictions of several Skyrme interactions in the RPA, with widths of 1 MeV introduced by 
hand following Ref. |24|. for the isoscalar-El strength distribution in 208 Pb. The figure also shows the locations of 
the measured low- and high-energy peaks of Davis et al. |2^| and of Clark et al. [l^- Nearly all self-consistent RPA 
calculations, including ours (except with SkP) over-predict the energy of the larger peak by a few MeV [2^, H^- SIII 
does a particularly poor job. The predicted low-energy strength is closer to experiment, though usually a little too 
low. Tabled summarizes the situation. Unfortunately, the data are not precise enough to extract much more than 
the centroids of the two peaks. Since the Schiff-strength distribution in 199 Hg helps determine the Schiff moment, it 
would clearly be useful to have better data, either in that nucleus or a nearby one such as 208 Pb. 

The isovector-El strength distribution also bears on the Schiff moment through the second term of the Schiff 
operator, but it is well understood experimentally and generally reproduced fairly well by Skyrme interactions. 

III. RESULTS AND ESTIMATE OF UNCERTAINTY 
A. Results with several forces 

The Schiff moment can be written as 



S = a gg + ai gg 1 + a 2 gg 2 , 



(12) 



6 





12 18 24 30 36 42 

Energy (MeV) 

FIG. 4: Isoscalar-El strength distributions for 208 Pb predicted by the Skyrme interactions SkP, SkO' and SIII in self-consistent 
HF+RPA. The experimental bounds on the low-energy (E x i) and high-energy (E X 2) peaks are also shown (see TabloQJ. 





low (MeV) 


high (MeV) 


SkM* 


11.0 


25.3 


SkP 


10.0 


23.4 


SIII 


11.6 


28.3 


SLy4 


11.4 


26.4 


SkO' 


10.3 


24.8 


Experiment [23J 




22.4 ±0.5 


Experiment [22J 


12.2 ±0.6 


19.9 ±0.8 



TABLE I: Comparison between experimental and theoretical results for the centroids of the low- and high-energy peaks in the 
distribution of isoscalar El strength in 209 Pb. The experimental results are from Refs. ETil and l23ll . (Ref. |23ll identifies only 
the high-energy peak.) The theoretical distributions are from self-consistent HF±RPA calculations with five different Skyrme 
interactions. 



where gi are the P- and T-odd 7rNN coupling constants and all the nuclear physics is summarized by the three 
coefficients a^. We present our results for these coefficients by showing the effects, in turn, of several improvements 
on early calculations. 

The first calculations of Schiff moments , as noted above, correspond to our first-order diagrams (iii) and (vi) of 
Fig. [2 but with no pairing, with a simple Wood-Saxon potential in place of a self-consistent mean field, and with the 
zero-range limit of the direct part of W . The results of Ref. are given here in the first line of Tab. |H] 

When we repeat the calculation, evaluating diagrams (iii) and (vi) with W approximated by its direct part in 
the zero-range limit and with the mean field from the Skyrme interaction SkO' (so that the only differences in the 
calculations are the one-body potential and BCS paring), we get the coefficients in second line of the table. 

The finite range of the potential reduces the a.; from these zero-range values by 30-40%, depending on the Skyrme 
interaction used. Exchange terms, when the range is finite, decrease ao by a few percent, have no effect on a\ and 
increase a 2 by half the amount they decrease ao. 

The three coefficients in lines 1 and 2 of the table are not independent; the isotensor coefficient is exactly two times 
larger than the isovector and isoscalar coefficients. Because the valence neutron must excite core protons to couple 
to the Schiff operator through (iii) and (vi) of Fig. ^ only the neutron-proton part of W contributes, and under the 
assumptions of no core spin, no exchange terms, and zero-range, W reduces to 

WKd*(r„-r p ) = -—4 (ggo + ggi + 2gg 2 )<r n ■ V n S(r n -r p ) , no core spin (13) 
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ao 


ai 


d2 


Ref. [6j 


0.087 


0.087 


0.174 


Naive limit 


0.095 


0.095 


0.190 


Diagram A only 


0.018 


0.034 


0.031 


Full result 


0.010 


0.074 


0.018 



TABLE II: Calculated coefficients ai from Ref. and with the Skyrme interaction SkO' in several limits (see text). The full 
result is in the last line. 

Thus, in the approach of Flambaum et al. , the Schiff moment is a function of a single parameter, usually called r) np . 
Exchange terms add another independent parameter and QRPA bubbles bring in a third. This last term arises because 
the valence neutron, besides exciting a proton quasiparticle pair in the core, can now excite neutron quasiparticle pairs 
that annihilate and create proton pairs inside the bubble that then couple to the Schiff operator (see Fig- EI) • Thus, 
in our complete calculation, particularly when the B diagrams are included, the three coefficients ai are independent. 

The collectivity of the core turns out to be very important. As can be seen in the third line of Tab. ITT1 when the 
single-particle bubble in the calculation above is replaced by the full QRPA bubble sum to give diagram A, all three 
<Zj shrink substantially. The reason is that the Schiff strength is pushed on average to higher energies, both in the 
low-lying and high-lying analogs of the isoscalar-El distribution. (The high-lying peak actually is replaced by two 
peaks, the higher of which is at about 38 MeV. There is no peak corresponding to the giant isovector-El resonance, as 
was shown in Ref. .29]].) The reduction is greater in the isoscalar and tensor channels — a factor of 4 to 6 depending 
on the Skyrme interaction — than in the isovector channel, where it is a factor of 2 to 3. Figure 03 shows the integral 
of the contribution to diagram A as a function of core-excitation energy (so that at large energy the lines approach 
the value for the diagram) for ai, with and without the bubble sum, for all 5 forces. 

The reason for the difference in the size of the reduction is that the go and g2 parts of the interaction affect 
protons and neutrons in opposite ways (see, e.g., Eq. (9) of Ref. ^3]), causing a destructive interference, while the 
gi part affects them in the same way. This difference is absent from the single-quasiparticle picture because neutron 
excitations of the core don't play a role there. Another way of saying the same thing is that when the neutrons and 
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FIG. 5: Integrated contribution of gi to diagram A as a function of core-excitation energy (the value of the isovector part 
of the diagram is given by the lines at large energy) in the two-quasiparticle limit (top) and in the QRPA (bottom) for five 
different Skyrme interactions. The two bumps are at the energies of the low- and high-lying peaks of the isoscalar-El strength 
distribution. 
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protons are affected in the same way, the second (dipolc-like) term in the Schiff operator, Eq. J2J), contributes very 
iittle because the center of mass and center of charge move together. In the other two channels the contribution of 
the second term is similar in magnitude and opposite in sign to that of the first term so that, as Fig.HOshows, the net 
value is smaller. 
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FIG. 6: Integrated contributions of the first (dashed) and second (dotted) terms in the Schiff operator, together with the total 
contribution (solid) to diagram A, as a function of core-excitation energy with the SkO' interaction. The top panel is for the 
isovector coefficient ai and the bottom for the isoscalar coefficient ao. 



The type-B diagrams (see Fig.[2J) are important corrections to diagram A. The effective weak isoscalar and isotensor 
one-body potentials (i.e. the tadpole) contribute with opposite sign from that of the isovector potential; again, see 
Eq. (9) of which displays the direct part of the one-body potential explicitly. The sign turns out to be opposite 
that of diagram A in the isoscalar and isotensor channels, further suppressing ao and 02, and the same as diagram 
A in the isovector channel, largely counteracting the suppression by collectivity in that diagram. The net result for 
SkO' is in the last line of Tab. [HI for the other forces the net results appear in Tab. IIIII The isovector coefficient a\ 
ends up not much different from the early estimate of Ref. |(| but the isoscalar coefficient ao is smaller by a factor of 
about 9 to 40 and the isotensor coefficient by a factor of about 7 to 16. 



B. Uncertainty and final result 



The several Skyrme interactions we use all give different results, but the spread in numbers is about a factor of four 
in the isoscalar channel, two in the isotensor channel, and much less for the large isovector coefficient a\. It is possible 
that all the interactions are systematically deficient, but we have no evidence for that. In any event, the effective 



ao ai a 2 

SkM* 0.009 0.070 0.022 

SkP 0.002 0.065 0.011 

SIII 0.010 0.057 0.025 

SLy4 0.003 0.090 0.013 

SkO' 0.010 0.074 0.018 



TABLE III: Full coefficients a^ in 199 Hg for the five different Skyrme interactions used here. The units are e fm 3 . 
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interaction is not the only source of uncertainty. We have evaluated only a subset of all diagrams, and although 
it is not obvious whether all the rest should be evaluated with effective interactions that are determined through 
Hartree-Fock- or RPA-bascd fits, there are some that should certainly be included and we would like to estimate 
their size. 

The diagrams labeled C and D in Fig. are the leading terms in bubble chains that would result from including W 
in the Hartree-Fock calculation (C) and in the QRPA calculation (D) in 198 Hg. We evaluated both sets in the simple 
nucleus 209 Pb (which has no pairing at the mean-field level), and found that diagram C can be nearly as large as the 
type-B diagrams in the isovector and isotensor channels. The same is true of diagram D in the isotensor channel. In 
that nucleus, however, diagram A is much larger than all the others and essentially determines the <Zj. In 199 Hg, we 
only evaluated diagram C, but found that even though diagrams A and B can cancel there, they do so the most in 
the isoscalar channel, where the diagrams C and D are smallest. In the end diagram C never amounts to more than 
10% of the sum of the A and B diagrams, and usually amounts to much less. Including the higher order (QRPA) 
terms in the bubble chain will only reduce the diagram-C contribution, so we conclude that it can be neglected. We 
are not positive that the same statement is true of diagram D, but unless it is much larger in 199 Hg than in 209 Pb 
(none of the other diagrams are), it can be neglected too. 

The diagram labeled E represents a correction from outside our framework that is of the same order as the terms 
we include. We evaluated it in 209 Pb; it is uniformly smaller than those of type C and D. Unless the situation is 
very different in 199 Hg, it can be neglected as well. The fact that these extra diagrams are all small is not terribly 
surprising; they all bring in extra energy denominators and/or interrupt the collective bubble. 
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(i) diag-C 



(ii) diag-D 



(iii) diag-E 



FIG. 7: Diagrams we did not include in our calculation in Hg but the value of which we estimated through calculations in 
209 Pb (and 199 Hg in the case of diagram C). We have omitted the labels on the the interactions; they are the same as in the 
earlier figures. 

We have also not included the normalization factor Af in Equation ©. When calculated to second order in 209 Pb, it 
is about 1.05, independent of the Skyrmc interaction used. Though this factor could be larger in RPA order because of 
low-lying phonons, most strength is pushed up by the RPA and we do not anticipate a large increase. It is reasonable 
to assume these statements are true in 199 Hg as well. 

At short distances the NN potential is strongly repulsive and the associated short-range correlations should be taken 
into account. Reference however, reports that the correlations reduce matrix elements of the effective one-body 
P- and T-violating pion-exchange potential only by about 5%, and in Ref. [3(j, which calculates the Schiff moment 
of 225 Ra, their effects are smaller than 10%. We are not missing much by neglecting them, though we would be if we 
included a p-meson exchange potential. 

When all is said and done, the uncertainty is dominated by our uncertainty in the effective interaction. Our 
preferred interaction is SkO', which (to repeat) gives the result 

Sg^'g - 0.010gg + 0.074g gl + 0.018gg 2 [e fm 3 ] . (14) 

If instead we average the results from the five interactions, we get 

SS c Hg - 0.007gg + 0.071g gl + 0.018gg 2 [e fm 3 ] , (15) 

The range of results in Tab. Illll is a measure of the uncertainty. 

As noted in the introduction, Refs. H, @ contain a similar calculation. They report 

SS^ 81 = 0.0004gg + 0.055g gl + 0.009gg 2 [e fm 3 ] , (16) 

the most striking aspect of which is the isoscalar coefficient a ; it is more than an order of magnitude smaller than 
our preferred value and five times smaller than the smallest coefficient produced by any of our interactions. We see no 
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fundamental reason for such serious suppression, and suspect that the same cancellation we observe is coincidcntally 
more precise in the single Hg calculation of Refs. H,ij. The authors applied their method to other nuclei, but did 
not find the same level of suppression in any of them. Even the cancellation produced in our calculations by SkP and 
SLy4 seems coincidentally severe. Though it is possible that other realistic Skyrme interactions would produce still 
smaller coefficients, we have a hard time imagining it. 

IV. CONCLUSIONS 

Our goal has been a good calculation of the dependence of the Schiff moment of 199 Hg, the quantity that determines 
the electric dipole moment of the corresponding atom, on three P- and T-violating 7rNN coupling constants. The 
current experimental limit on the dipole moment of the 199 Hg atom, \d\ < 2.1 x 10 _28 e cm, together with the 
theoretical results of Ref. d = -2.8 x 10 _17 (S'/e fm 3 ) e cm, yields the constraint \S\ < 7.5 x 10~ 12 e fm 3 . The a, 
calculated in this paper then give a constraint on the three gj. 

In obtaining the a t we have included what we believe to be most of the important physics, including a pion-exchange 
P- and T-violating interaction, collective effects that are known to renormalizc strength distributions of Schiff-like 
operators, pairing at the mean field level, self-consistency, and finally, several different Skyrme interactions. The last 
of these, together with an examination of effects we omitted, allows us to give the first real discussion of uncertainty 
for a calculation in this experimentally important nucleus. 

We conclude that while the isovector coefficient a\ is not very different from the initial estimate of Ref. @, the 
isoscalar coefficient, which determines the limit one can set on the QCD parameter 6, is smaller by between about 9 
and 40 (with the former our preferred value) and the isotensor parameter ci2 by a factor between about 7 and 16 (with 
our preferred value about 10). The uncertainty in these numbers comes primarily from our lack of knowledge about 
the effective interaction. There is good reason to make better measurements of low-lying dipole strength, particularly 
in the isoscalar channel. They would help to unravel the details of nuclear structure that determine the Schiff moment. 
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V. APPENDIX 
A. Rules for quasiparticle diagrams in the uncoupled basis 

There are some differences between our rules for quasiparticle diagrams and the usual rules for particle-hole dia- 
grams. The main one is that one- and two-body operators are written in a quasiparticle basis and do not conserve 
quasiparticle number, leading to different expressions for matrix elements. An example is the generic quasiparticle 
operator O20, which contains two quasiparticle creation and no destruction operators. Its matrix elements will be 
written (fcZ | O20 1 — )> which means that it creates two quasiparticle states \k)\l) out of the quasiparticle vacuum |— ). 

In what follows, "in" refers to lines with arrows pointing toward the vertex and "out" to lines pointing away from 
it. A diagram should be read from top to bottom, and from left to right. The rules are then: 

1. Each operator On contributes (out | On |in); 

2. Each operator O20 contributes (out, out'|02o| — )j because the diagram is read from the left, the label "out" is 
on the line that is further to the left; 

3. Each operator O02 contributes (— |Oo2lin,in')i 

4. Each operator O22 contributes (out, out'|022|in, in'); 

5. Each operator O31 contributes (out,out',out"|03i|in); 

6. Each operator O13 contributes (out|Oi3|in, in', in"); 
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7. Each operator O40 contributes (out,out',out",out"'|04o|— ); 

8. Each operator O04 contributes (— |Oo4 |in, in', in", in"'); 

9. The diagram should be summed over all intermediate states; 

10. Energy denominators are evaluated by operating with (e a — Ho) between the action of every two operators 
in the diagram, giving [e a + J2 k e k] ■> where e k are quasiparticle energies; 

11. The phase for each diagram is (—)"', where ni is the number of closed loops. 

12. A factor of 1/2 is included for each pair of lines that start at the same vertex and end at the same vertex. 
Folded diagrams with additional rules occur in general, but we do not discuss them here. 



B. Matrix elements of quasiparticles operators 

We first summarize some important quantities involving the quasiparticle creation and annihilation operators 
and q, which are defined in terms of the usual particle operators a' and a by 

f q k = u k a k - v k a\ 
1 q\ = u k a\ - v k a k 



(17) 



<lk = u k a k + v k a\ 

Ok = U k^k + v kO-k 



Here 



qk = qi kjkmk = H lk+lk+mk qi kJk ~m k = {-) l ^» +m *q- k . (18) 

From the anti-commutation rules for a* and a, we derive the following anti-commutation rules for the quasiparticle 
operators in Eq. (|17|l 

Uk,qi} = ULqj} = {qk,qi} = {qi,qj} = {qk,qi} = {q\,q\} = 

{qk,q\} = {qkA\} = &kl (19) 
{q k ,ql} = (-Y k+Jk+mk 5-kl 

Using definition i|17f) and properties l|18|) and l|19|l . one can write a one-body operator in second quantization as 

T = Tu a\ai =T +T n +T 20 + T 02 , (20) 

kl 

where 

n = J2 v l T ^ ■ ( 21 ) 

k 

Tn = ^T k i(u k uiqlqi - v k viqjq k ) , (22) 

kl 

T20 = y^T k iUhViq\.qi , (23) 

kl 

T 02 = ^T k iuiv k q k qi . (24) 



kl 



Here, the subscripts indicate the number of quasiparticle creation and annihilation operators involved. In the same 
way, a two-body operator takes the form 



klij klij 

= V + Vu + V 2Q + V 02 + V22 + Vm + Via + V 40 + V04 , (25) 
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where 



= \^{ukVkumPkPiV k -ki-i+2vlv?Vkiki) , (26) 

kl 

Vu = - ^2 ukViUiViPiVku-i (q\qi + qjqk 



2 

kli 



+ ^2 v k v kiki (uiu i q\q l - viViq\qij , (27) 
kli 

7 ^2 u i V i P i V kli-i ("fe"/?!^ + *W2M) 



^20 = 

+ ^ Va^M, (28) 

kli 
klij 

+ ^2u k viU j v i Vkii J qiqUi ( lj > ( 29 ) 

klij 



V 31 = ^^UkViVkiij (uiUjqjqjqjqk + ViVjqlqtqlqi^ , (30) 

klij 

V 40 = j^UkUiViVjVkUjqlqjqlql , (31) 



klij 

with P k = (_)l»+i*+"»* J = (y ot )t and V abcd = V abcd - 

The matrix elements of quasiparticle operators are related to those of the usual one- and two-body operators. We 
show how this works for the operator T20; the generalization to other operators follows automatically. T20 adds two 
quasiparticles to any state. From \ij) = q\q\\— ) and Eq. lt25|) . we write 

(ij\Tm\-) = ^T k iu k vi(-\q jqi qlq\\-) . (32) 
kl 

Since 

(-\QjqM\-) = Pti-SjkS-u + S^Sik] , (33) 

Eq. lO becomes 

(*i|Tao|-> = -P l u J v t ( 3 \T\-i)+P J u l v J (i\T\-j). (34) 
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